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Abstract . Compositional variations in nickel , iron, and cobalt in nine octahedrites 
were determined by electron microprobe X-ray analysis . The steep gradients in 
taenite extending up to 35-50 wt.% Ni and the shallow decreases (0.2-1 wt.% Ni) 

„ in kamacite near the interface of the kamacite and raenite phases were formed 

during the final stage cf formation of the Widmanstatten structure as cooling prevented 
attainment of equilibrium by diffusion . A simplified theoretical model based on 
the Fe-Ni system was used to relate the extent of diffusion to cooling rates of the 
meteorite parent-bodies . These depend critically on assumptions concerning the 
nucleation temperatures of kamacite, pressure, phase equilibria, and cooling model . 
Metallurgical evidence implies that nucleation of kamacite may not have occured 
at the y / ct + y boundry in the Fe-Ni system, but after undercooling of about 
I00°C at I atmosphere pressure. d z nuclei may have initiated diffusion growth 
of kamacite which ceased by 350°C . Differing nucleation temperatures account 
for meteorites of different structure, but identical nickel content, for a range of 
! kamacite band-widths, and for coarse plessite . Fine plessite, of variable composition 

(10-28 wt.% Ni) developed after partial transformation of y to d 3 below . It 

seems most reasonable that the Widmanstatten structure formed at low pressures in 

7 8 "6 

10 -10 years at cooling rotes of 2-l2°C/l0 years in the temperature interval 650- 

350°C . These are consistent with a chondritic parent-body (or bodies) with a metallic 

core and a radius of l00-200km . 




INTRODUCTION 


The existence of high-nickel borders in the taenite phase of iron meteorites 
has been shown by previous \electron microprobe X-ray analyses (Yavnel 1 , et. al., 

1958; Maringer, et. al., 1959; F eller-Kniepmeyer and Uhlig / 1961). These non- 
equilibrium nickel gradients are related to the formation of the Widmanstatten 
structure by phase precipitation and growth by solid state diffusion . Uhlig (1954) 
and Urey (1956) have tried to estimate the time for formation of the Widmanstatten 
pattern by using the simpliest diffusion parameters, as did Massalski and Park (1962) 
by use of the Fe-Ni phase diagram. The nickel diffusion gradients in the taenite 
phase of nine octahedrites, have been measured with an electron microprobe . A 
refined diffusion theory and cooling model was applied in order to determine the 
probable cooling rates and sizes of the meteorite parent-bodies . 

EXPERIMENTAL METHODS 

The meteorites were selected so as to provide a *ange of kamacite band- 
widths in all octahedrite classes and to represent all four Ge-Ga groups ( Lovering , et .al ., 
1957). Octahedrites of different texture but similar nickel content were chosen — 
e.g., Bristol and Toluca — to test the hypothesis ( Goldberg, et. al., 1951) that 
different band-widths might be due to different cooling rates . 

Unweathered sections were cut with a carborundum wheel, polished with 
silicon carbide and diamond paste through I p, then cleaned ultrasonicaliy in 
benzene. No etching was done, in order to avoid surface redeposition (Feller-Kniepmeyer , 
1961). The taenite bands and plessitic areas to be analyzed were selected by 
microscopic examination. Generally, the larger symmetrical bands of taenite or 
plessite between the primary kamacite bands defining the Widmanstatten structure 
were chosen . Particular care was used to avoid areas near intersecting band systems 
where diffusion had proceeded from two or more directions . Irregular areas or bands 
near inclusions were avoided. Figure i illustrates a typical band. 

Analyses were made at the Scripps Institution of Oceanography, with a 
modified electron microprobe X-ray analyzer (Applied Research Laboratories, Glendale, 
California). The instrument was operated at 25 kv. with a probe diameter of less 
than I jj . Appropriate integration times and stop-counting gave precisions in measure- 
ments of 1 0.2%. The microprobe tracks were made perpendicular to the long surface 
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axis of the band near the center, as shown in Figure I . The probe was stepped 
1-9 p between analyses., depending on the compositional gradient. Microscopic 
calibration confirmed the accuracy (0.2 p) o f the microprobe stepping mechanism . 

Pure metals were used as standards for iron, nickel, and cobalt. A 
chemically analyzed schreibersite (15.5% P) in the meteorite Horse Creek was the 
phosphorus standard . All raw date were corrected foi detector dead-time and back- 
ground . Mass -absorption corrections were applied to nickel, iron, and cobalt; 
fluorescence corrections were applied to iron . Secondary fluorescence of cobalt by 
nickel K -j^ is negligible . The accuracy of the theoretical corrections was checked 
by comparison with chemically analyzed Fe-Ni alloys (International Nickel Company; 
furnished by J. Goldstein, Department of Metallurgy , Massachusetts Institute of 
Technology) . The results are given in Table I . 

In the alloys the microprobe accuracy for nickel was better than± 0.5%, 
and the accuracy for iron± 1%. In meteorite areas where t.ie composition changed 
across the diameter of the X-ray producing volume, the accuracy was less . In fine 
plessite, unresolvable by the microprobe, only average compositions could be 
determined, but in tcenite, where the composition varied smoothly, an accurate 
concentration profile was easily obtained. In the areas of the meteorites where 
rhabdite o r schreibersite was not present, the measured phosphorus content was very 
close to the detection limit of 200 ppm. 

The orientation of the microprobe track relative to the [[lo) ^ zone axis 
was determined by taking Laue photographs (molybdenum radiation) of the adjacent 
kamacite . Mehl and Derge (1937) have shown that during formation of the 
Widmanstatten structure the oC-(IIO) planes of kamacite precipitate as plates parallel 
to the (1 10) pianes of V i.e., in the (X-phase the direction of crystal growth by diffusion 
is the [llqj* zone axis . The angle between this direction of diffusion and the micro- 
probe track was determined by indexing the Laue photographs and determining the 
corresponding angle between the direction and the specimen nonna! . 

RESULTS 

Nickel profiles . Typical profiles of the nickel concentration across selected taenite 
or plessite bands are shown in Figures 2-11 . The scales permit comparison between 
meteorites of the steepness of the concentration gradients, the distance between 
the phase boundry and the nucleation point (see Appendix) the smoothness of the 
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gradient, and the central and bulk nickel concentration . 

Most bands found in coarse and medium octahedrites (Figures 2-7) and narrow 
taenite bands found in fine octahedrites show the interior nickel content is higher 
than the bulk nickel content because nickel has diffused from the boundry to the 
center. Wide plessitic areas in fine octahedrites have the same average central 
nickel contents as the bulk concentration for the meteorite (Figures 8-10). 

Symmetrica! gradients across the entire taenite band were found only (e.g.. 

Figures 2 and 6) when equal widths existed in the adjacent kamacite bands and the 
widths of the outer taenite bands were almost the same . When the above condition 
was not present profiles were usually observed, with the minimum nickel concentration 
in taenite lying off-center toward the narrower kamacite band (e.g.. Figures 4,5, 
and 7) . Assymmetric profiles were also observed when tracks were made at large 
angles relative to the (Tioj direction of kamacite or were near intersecting band 
systems such that diffusion had proceeded from more than one source . The central 
nickel content of taenite bands within one meteorite varies with the taenite band- 
widths, narrow bands showing proportional!^ higher nickel (e. g.. Figures 4 and 5) - 
Wider kamacite bands introduce more nickel into the bordering taenite, flattening 
the gradients at the border and raising the central taenite content (Figures 4 and 7) . 

Figure 7 demonstrates the effect on the nickel -profiles of varying kamacite and taenite 
band-widths . The profiles are least affected near the boundry; in fact, they 
merge into a common gradient as they approach the boundry. 

The kamacite was remarkably homogeneous except within 1 0-30 of the 
<X/Y border . A slight structural dependence of mean nickel content of kamacite was 
observed: 6.6 wt. % for coarse octahedrites, 7.0 wt . % for medium octahedrites 
and 7.3 wt . % for fine octahedrites . At distances less than 10-30 jj from the boundry 
a continuous slight decrease in nickel was observed without exception as the boundry 
was approached, in agreement <vith observations by YavneTet .a I .(1958) and Agrell ; et .a l ., 
(1963). The decreases ranged from about 1% for coarse octahedrites (Figures 2-4) 
less for medium octahedrites (Figures 5-7) to only 0.2% - 0.3% for fine octahedrites 
(Figures 8-10) . The extent of the decreases was correspondingly less . The decrease 
is real and not due to mass-absorption effects or sample preparation because: (a) it is 
independent of the direction of movement of the probe; (b) it is present whether or 
not the surface was finely or roughly polished; (c) it is present at all of the many inter- 
faces analyzed; (d) taenite bands with steeper edge -gradients showed correspondingly 
shorter depression in kamacite . 
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Plessite . Microprobe analysis of nlessite showed the same compositional variation 
that Perry (1944) observed. Coarse plessite, such as perthitic lamellas in coarse 
octahedrites, was resolvable with the microprobe into distinct taenite and kamacite 
regions; the fine plessite often occuring in fine octahedrites (Figures 8 and 9) was 
not resolvable; plessite of intermediate size (Figure 6) occurs frequently in medium 
octahedrites . Moving from the boundry to*vard the interior of taenite, the smooth 
nickel gradient starts showing fine structure of progressivily greater compositional 
variation at lower nickel concentrations (Figure 8 and 12), while the average plessite 
content follows the diffusion gradient until the interior is reached. Plessite in fine 
octahodrites had the same average composition as the bulk chemical content, pro- 
vided enough microprobe analyses were made in plessite to give good statistics . 
Meteorites of fine structure have finer plessite than coarser octahedrites and have 
smaller amplitudes of nickel variation in plessite, as seen by comparison of Figures 
6 and 8. The average nickel content of the kamacite in fin<_ plessite approached 
about 6 wt. % and is lower than the nickel content of the kamacite in coarse plessite 
or the Widmanstatten structure . 

The &/i boundries in plessite, as seen in Figures 12 and 13, were not in 
straight octahedral orientation as was the case in the large Widmanstatten structure . 
Rather, warped boundries with preferred orientation were comm'' . The existence 
of plessite does not always correlate simply with nickel content, ,ven in one meteorite . 
Figure 3 shows a taenite band of which only one side has exsolved into plessite, while 
Figure 4 shows two taenite bands, of which the one with the higher nickel content 
has exsolved. The upper limit in nickel content for plessite occurence varies widely 
from approximately 14-28 wt . % nickel . There appears to be no lower limit for 
plessite formation; Bristol (Figure 10) has only 8.3 wt . % nickel, but has many areos 
of fine plessite . 

Other elements . Iron was determined concurrently with nickel and invariably varied 
in an inverse manner, such that Fe + Ni S 100% . Cobalt was measured in many bands 
and followed iron, with a «/V distribution coefficient favoring kamacite, in agreement 
with observations by Nichipo ruk (l 996) .Figure II shows a track across a taenite band in 
Canyon Diablo in which a pair of kamacite lamalle have apparently exsolved in the 
late stages of formation of the nickel gradients . The inverse correlation between 
nickel and cobalt is evident. 

The phosphorus content of the metal-phase was at the detection limit (0.02%) 
in all meteorites except Canyon Diablo, in which infrequent rhabdite inclusions 
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occured and which had a central phosphorus content in kamacite 0.1%, but which 
decreased as taenite bands were approached. 

Orientation studies . Table 2 lists the angle 90-# between the track direction and the 
growth direction . A correction factor c os O was applied to the measured profile 
distances (x 1 ) for diffusion calculations . In three cases (marked by *) subgrain- 
boundries had apparently caused imperfections in the crystal structure, the Laue 
photographs were diffuse, and no orientation could be determined. Such an effect 
is apparently independent of the perfection of the Widmanstatten structure (Derge 
and Kommei , 1937). 

DISCUSSION 

Phase transformations in the Fe-Ni system . The classic and generally accepted ex- 
planation for the formation of the Widmanstatten structure was offered by Derge and 
Kommei (1937) and Perry (1944). A homogeneous, solidified f c.c. phase i cooled 
slowly until the low-nickel b.c.c. phase (oc or kamacite) nucleated and grew by 
solid state diffusion, absorbing iron and expelling excess nickel into the parent-phase 
taenite. The (110) crystallographic plane of the ^ -phase precipitated parallel to the 
(III) octahedral planes of the original V-phase ( Young , 1926; Mehl and Derge, 1937). 
Low-nickel hexahedrites had all of their Y -phase absorbed at lower temperatures as 
the meteorites bulk concentration passed back into the one phase region oM stability. 
The Fe-Ni phase diagram (Figure 14) approximates this sequence of events in general . 
The slope of the K/a^+Y boundry explains the general decreasing width of the 
Widmanstatten structure seen in Table 3 as the bulk nickel content increases . As 
the kamacite grew into the taenite, the interface concentrations at the moving 
boundry were dictated by the equilibrium phase compositions. The inhomogeneties 
seen in Figures 2-11 (gross in taenite; fine in kamacite) show that equilibrium was 
not maintained away from the kamacite -taenite interface. Rather, diffusion slowed 
down and was finally arrested at lower temperatures, leaving the gradients in nickel, 
iron and cobalt as evidence of a cooling p/ocess . These compos i Mona I gradients are 
quantitatively related to the cooling rate and times of formation of the Widmanstatten 
structure, which may be determined, provided the correct phase equilibria, nucleation 
and diffusion theory, and cooling model for the parent-body are assumed. The diffusion 
coefficient of nickel inj^, the total nickel content, the geometry and temperatures 
of nucleation of the cooling rates of the meteorite at these temperatures and 
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pressure as it effects these variables are the parameters which determine the size 
and form of nickel -profiles in the Widmanstatten structure. 

The currently accepted Fe-Ni sub-solidus phase diagram (Figure 14) of Owen 
and Liu (1949) explains the formation of the Widmanstatten structure in general, al- 
though the nickel content of kamacite in octahedrites is higher than predicted. This 
may be because th e + V boundry lies to the right or because the boundry is shifted 
by the presence of a few tenths percent of phosphorus in solution (Ringwood and Kaufman , 
1962). The 0.5% cobalt content will not shift the phase boundries significantly. Owen 
and Liu's (1949 )/ /<*»/ boundry and Goldstein's (1964) ientative «/<* +• boundry uie 
assumed to approximate the octahedrite systems . The decreases in the nickel content of 
kamacite near kamacite-taenite boundries are most easily explained if the <*/ + v 
boundry in the meteoritic phase system moves to less nickel content at lower temperature, 
(400-350°C) where boundry movement ceases and kamacite cannot homogenize . 

Vogel (1925, 1951, 1937 ) has hypothesized that the Widmanstatten structure 
formed during a high temperature exsolution of <*, as dictated by the Fe-Ni -P system. 
However, this seems unlikely, since the phosphorus content of octahedrites is too low 
(0-0 .5% P); also no correlation of structure with phosphorus content exists . Even more 
significantly, it is calculated by Equations I and 2 that even as low as !IOO°C, homo- 
geneous _/-phase of 5 cm can be formed . 

The kinetics of sub-solidus phase transformation in the Fe-Ni system have been 
misunderstood f.equently in the past. It is of critical importance in interpretating the 
Widmanstatten structure to define the nucleation temperatures of kamacite. This is 
probably not the temperature at which the bulk composition ciosses the /A* Y 
boundry, but at some lower temperature . Cooling experiments on crtificial alloys by 
Jones and Pumphrey (1949), Allen and Earley (1950) arid Kaufman and Cohen (1956), 
as well as earlier heating experiments cf Hanson and Hanson (1920) and Owen and Su ll y 
(1939) have shown that, upon cooling, _/-phase undercools below the //«<•*• v' phase 
boundry. At the temperature the f.r.c. V'-struc'-ire starts to change by a diffusion- 
less transformation to^, a strained b.c.c. structure of the same composition. Anrooitrig 
or very slow cooling permits <*2 f° decompose into in accordance with the phase 

diagram . 

Kaufman and Cohen (1958, 1956) have reviewed the thermodynamics of tht: 
transformation and showed that at T q AF(V -^e^) = 0. Further undercooling tn Fe-Ni 
alloys from T q to is required because of interfacial and strain energy contributions 
to the free energy. The transformation is not complete until a still lower temperature 
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My is reached . The and T q temperatures shown in Figure 14 are not accurately known 
and may vary. m ay be raised or lowered by impurities; cobalt increases M s lightly 
(Yeo , 1963). Stress increase s M & b y IO-l5°C/kb., and pressure lowers by 8°C/kb . 

Patel and Cohen ,1953 ). Likewise, pressure lowers the equilibrium phase boundries 
(Fling wood and Kaufman , 1962 ), and phosphorus raises them ( Vogel, 1951). However, for 
reasons given by Ringwood and Kaufman (1962)and the diffusion calculations below, 
minor components in meteorites do not seem to have affected the phase boundrie» sign! — 
ticcntly . Thus the Y// + * boundry Qnd i n the Fe-Ni system may be taken as 

giving approximate upper and lower limits on the temperature at which bc.c. formed in 
iron meteorites . 

If one extrapolates the studies on artificial alloys to geologic time intervals (an 
assumption that does not contradict present experimental evidence), undercooling of 
as far as is possible . Given sufficient time at temperatures oetween T^and M e , <*. 


nuclei might exceed the activation barrier for formation and quickly grow as a very thin 
(p-size) plate parallel to the (lll)^plane, stopping laterally only when a grain boundry 
or another <^2 was encountered. This b.c.c. phase would then be the nucleus for de- 
composition into the equi librium or V . The assumption of different degrees of under- 
cooling is supported by the data in Figures 2-11, since kamacito bands of widely vary- 
ing width are seen in the same meteorite; these have expelled differing amounts of nickel 
into taen'te. The assymetric taenite bands with different nickel -distance profiles are also 
explained bv different nucleation temperatures for the adjacent karpacite bands . Local 
compositional differences are probably too small to account for these effects, since the 
average microprobe analysis of fine plessite (e.g., Altoncih) is the same as the bulk compo- 
sition. Meteorites of different mean kamacite band-width, but identical nickel content 
(see Massalski , 1962) could also have had different average nucleation temperatures, 
rather than different cooling rates, os Goldberg ,et .al . (1951) hypothesized. 

Since the gross Widmanstatten pattern suggests nucleation and grov> h of _f*. 
in / rather than growth of V in the assumption most consistent with studies on 
artificial alloys is that nucleation occured some place between T^ and M , the temp- 
erature at which becomes thermodynamically stable . Once formed, the ^^ 
nuclei would permit normal V -phase segregation . Bowles and Barrett (1952) have 
suggested that the martensitic mechanism may be operative in the early stages of 
nucleation and growth from solid soluiion . The mechanism proposed here differs 
somewhat from Owen's (1940) hypothesis that the Widmanstatten structure formed by 
decomposition ofo^below in that * „ nucleates the growth of<X in / . Random 
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nucieation at different temperatures between T q and would then permit kamacite 
bands of different widths to develop in the same meteorites or in meteorites with the 
same bulk composition but different structure . Nucieation of Widmanstatten growth 
at M seems unlikely , since all ^-phase (<30% Ni) would have instantaneously 


transformed to o^, requiring 1 to grow in**^> 


An explanation for the genesis of the Widmanstatten structure is incomplete 
without accounting for the origin and nature of plessite . The coarse p!es?ite consisting 
of narrow lamellas of kamacite and taenite is simply understood in terms of lower 
temperature nucieation off* from (Derge and Kommel / Iy37) . This occurs naturally 
in bands of taenite with relatively iow nickel content (e.g., Figure 6). Likewise the 
octahedral ly oriented needles of kamacite which are often observed in the cei tar of 
larger taenite bands are also understood as having formed like the large Widmanstatten 
pattern , but at lower temperatures with slower diffusion rates because of the higher 
interior nickel content of tlie taenite - Furthermore , the presence of kamacite spindles 
with octahedral orientation in nickel-rich ataxites (Ferry , 1944) shows that the mechanism 
of Widmanstatten -pattern formation extends into the plessitic nickel-rich ataxites . 

This implies that Uhlig^ (1954) cut-off of the Widmanstatten structure at 13% nicUel 
(17%, Table 3) Is meaningful only as a large-scale textural description, and the same 
mechanism operates at higher nickel contents and correspondingly lower temperatures, 
but on a necessarily smaller scale . The narrow kamacite plates does not extend 
throughout nickel -rich ataxites because the nucieation was localized. 

The formation of fine plessite seems to require a different mechanism than 
merely lower nucieation temperatures due to higher nickel content. The phase orienta- 
tions seen in the fine plessite of Figures 12 and 13 lack sharp parallel octahedral 
boundries. The random micro-taenite particles shown in Figures 8 and 9 suggest the 
possibility that they grew from kamacite, rather than the reverse. The dimensions of 
plessitic phases are discontinuous with the general dimensions of the Widmanstatten 
structure; this would nor be expected in fine octahedritei (Figures 8-10) if a range of 
nucieation temperatures were possible. It is suggested that Owen’s (1940) hypothesis 
of the origin of plessite is applicable to the formation of fine plessite . Undercooled 
V -phase evenUnally reached the M temperature (Figure 14) and started to be trans- 
formed to b.c.c.,^ 2 'Vb° $e > which then decomposed into + ^ , provided the temp- 
eratures were high enough to permit diffusion . The temperatures would be too low 
when transformation took plac* in the high nickel borders of taenite . Imperfect 
phase orientations would be expected. The mechanism is in agieement with laboratory 
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studies of the kinetics of Fe-Ni alloys and with the widespread occurence of fine 
plessite in meteorites . It also explains the occurence of nickel -rich ataxites of 
low nickel content. The general textural uniformity of fine plessite within fine 
octahedrites also suggests widespread nucleation ot a particular temperature dictated 

by 

Diffusion theory . The solutions of the diffusion equations to represent the concentra- 
tion gradients shown in Figures 2-10 :5 a difficult problem, involving non -isothermal 
phase growth with changing boundry conditions . it has not been attempted to derive 
an exact solution for the concentration-distance profiles, but rather, to derive 
appropriate diffusion parameters which may be related to thecooling rates and times 
of formation of the Widmanstatten pattern and subsequently to parent-body sizes of 
the meteorites upon assumption of an appropriate cooling model . The kamacite band- 
widths and the amount of nickel transferred across the <*/ * interface during phase 
growth are critically dependent upon the temperature of nucleation, the diffusion 
coefficient and the cooling rate of the meteorite within the temperature range of 
boundry movement. Because nickel inV is the slowest diffusing species, ( Hirano , 
eh al v 1961) D (nickel, Y ) controls the rate of &/Y boundry rrv vement. Recent 
accurate ir.terdiffusion measurements of D (nickel, Y ) by Goldstein et.al , (l°64), which 
may be extrapolated to 500°C with only + 50% error, were used to calculate tne 
cooling rate r during Widmanstatten pattern formation by use of the expression 
(Armstrong , 1958) 


D RT 2 
'v n n 
r Z : 

F E 


( 1 ) 


where D n is D at the nucleation temperature T^, r is the linear cooling rate, E is the 
activation energy for diffusion, and_F is a diffusion parameter (see Appendix) calcu- 
lated from 
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( 2 ) 



5 Is the amount of kamacite boundry movement sincej* -nucleation, and .T is evaluated 
from the boundry conditions . 

The possible cooling rates for each meteorite were calculated at four temperatures, 
including M > , T q , and the temperature at which the bulk concentration intersects 
the *7°** Y boundry . These are plotted against l/T in Figure 15 . 

The general slope is given by the activation energy for J) (nickel ); minor 
variations are due to differing C^and f . The meteorites are listed in order of de- 
creasing coarseness of structure, and different j5and_F values account for the spread 

of the lines . The average cooling rate for nucleation at T for seven octahedrites 

o 4 ~ o 

(excluding Bristol and Altonah) is 0.04 C/10 yrs . Bristol and Altonah either cooled 
somewhat faster or undercooled closer toM^. The identity of their cooling rates, 
diffusion gradients (compare Figures 9 and SO), and chemical composition ( Goldberg , 
et. al. , 1951; Nichiporuk ond Brown , 1964) suggest that they are genetically related. 

The particular band or bands studied in each meteorite have unique nucleation 
temperatures and band-widths, which are not necessarily representative of the mean 
(Table 3) band in the meteorite, nor in its class (Og, Om, or Of). The mean nickel 
contents and kamacite band-widths of octahedrite classes (Table 3) were obtained by 
calculating averages of the "superior" modern r.ickel analyses by Goldberg, et. al., (1951), 
Lovering et al ( 57)and others and of kamacite band-width measurements of Lovering (1956) 
and Goldberg, et. al ., (1951). t is given by one half the mean band-width corrected 
by a factor of 0 .8 to compensate for band width measurements made on meteorite 
sections which were not paralle' the original direction of diffusion . The correspond- 
ing cooling rates for coarse, medium and fine octahedrites are also plotted in Figure 15. 

At T^ they are almost the same, averaging 0 .12°C/I0^ yrs . Assuming that 
the degree of undercooling did not vary widely for different meteorites, it seems likely 
that the cooling rates during Widmanstatten -pattern formation are very nearly the 
same for all octahedrites . This condition implies that octahedrites were once at the 
same depth in a single parent-body or had the same cooling environment at 600°C 
in several bodies . 

Octahedrite parent-bodies , The nature of the meteorite parent-bodies has been 
discussed by many workers ( Anders , 1961 , 1964; Fish, et. al., I960; Goles, et. al. , (960; 
Lovering, 1957 a, b, 1962; Mason, I960, 1962; Ring wood , 1961; Urey , 1958, 1959, etc.). 

One of the simplest cooling models assumes tliat iron meteorites cooled from the melt 

in the center of a differentiated planet with a chondritic mantle . The relation (see Appendix) 

between the central cooling rate and the radius of the body is shown in Figure 16. 
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Long-lived radioactivity (K^®, U^, Th^'^) deci\ ises cooling rotes somewhat 

in bodies larger than 200 km . Figure 17 shows the corresponding parent-body radii for 
different nucleation temperatures; the comments made on Figure 15 apply also to figure 
17. For nucleation at T , the average radius for seven octahedrites, is 160 km and 
Sg, Om, and Of octahedrites could have cooled in a body of 100 km radius . For 
nucleation between T q and a 170 krn body would account for all twelve cases con- 
sidered, For example, an average medium octahedrite, cooling at the core of a chondritic 
body of this size at a rate of 2 .5°C/!0^ yrs . , would undercool !30°C to 590°C before 

nucleating, and formation of the Widmanstatten structure would be complete in less than 

8 8 
10 yrs . The chondritic body would have required about 3 J 10° yrs . to cool from I800°K 

to 350°C . 


Pressure is a parameter which must be considered, because of its effect on phase 

boundries and diffusion rates . High pressures over 10 kb. significantly depress the //a f * 

boundry in the Fe-Ni system ( Ring wood and Kaufman / 1961). Uhlig (1954) used this effect 

to explain the rarity of the Widmanstatten structure >13% nickel . However, the existence 

of undercooling of about !00°C and the occurence of micro-Widmanstatten patterns in 

nickel-rich ataxites (Perry, 1944) make this assumption of high static pressure unnecessary . 

It can be shown by use of Equation 0) (Appendix) that planetary-bodies greater than about 

400 km radius (with central pressures of 4 kb .) cannot even have cooled from l800°Kto 350°C 

9 

at the center in 4.5 4 10 yrs . This was also pointed out by Allan and Jacobs (1956) and 
Urey (1956). Furthermore, Goldstein (1964) has shown that static pressure strongly de- 
creases the diffusion coefficient of nickel in jjj-phase, such that it appears impossible to 
form coarse Widmanstatten structures under high pressure, even with no undercooling. 

The possibility can also be considered that the cooling rates calculated here 
could be due to octahedrites cooling in pockets near the surface of larger bodies, as en- 
visioned by Henderson and Perry (1958) and Urey (1959) . Even so, low pressures are required . 
Anders and Goles (1961) state other arguments against such a model . 

Perhaps more complex combinations of pressure, phase equilibria, and cooling 
model might be devised as our understanding of iron-pallasite-chondrite relationships 
grows . For example. C arter and Kennedy (1964) suggest that parent-bodies capable 
of producing diamond might have disrupted at temperatures above 800°C , then cooled 
more rapidly as the Widmanstatten structure developed. However, this hypothesis seems 
inconsistent with the strong general correlation of nickel -content with kamacite band- 
widths, and the similar cooling ratejcalcuiated for the meteorites studied. 
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APPENDIX 

Diffusion theory . Solutions of the diffusion equation for a constant diffusion coefficient 
D may be used for time-dependent diffusion coefficients by employing the parameter 
(Crank, 1956) defined by 


F = 



D( t ) dt 


0 ) 


For the case of a linear cooling rate r, Armstrong (1958) derived an approximate expression 
for F 


F = 


D n RT n 


rE 


( 2 ) 


D n is D at the nucleation temperature r i s the cooling rate, and_E_is the activation 

energy for diffusion . By combining Equation (2)with expressions for T (t) and D (f), one 

can derive a relation for the time t * for formation of the Widmanstatten structure between 

n-r 

and Tj, a temperature at which lower diffusion rates cause negligible boundry movement 


f n-f 


F T f 

D n T n 


Dn _ F E(T n -T f ) 
Df D n R T n 2 


In 


( 3 ) 
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Jost (1952) has given analytical solutions to *he diffusion equation for isothermal 
moving-boundry problems in two phases. For the purposes of approximation, we will 
assume by analogy with dost , that the following relation is sufficiently correct for linear 
cooling rates over fairly narrow temperatures ranges: 


C = 2 X /f” 


( 4 ) 


is the amount of boundry movement andj[ is a dimensionless parameter dependent upon 
the boundry conditions. Since the observed kamacite is almost homogeneous, .£ may be 
determined by 


£'«V c d) = fl 

C - £c0S 9 


where A is the total amount of nickel expelled into the V-phase . is determined by 
integration of the nickel -concentration-distance profile in taenite abo/e to the 
minimum in the profile . If the boundry conditions are constant,^. may be evaluated 
Uost, 1952) from the expression 


= fW Y e 


Q 4- ^ 


erfc y 


( 6 ) 
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where Q= (Cj| |). C|| and C| are the respective nickel concentrations at 

the interface between the V and ot-phases . Figure 18 shows that Q can vary greatly 

between T and T £ , but V varies much less . Because of the exponential decrease of 
_JL — L — 

JD with temperature., most of the phase growth and material transport across the inter- 
face occurs at and just below T n , with relatively small changes in boundry conditions; 
Therefore, Equation (4) may be assumed to hold approximately, using a mean Y , Y . 
Assuming that one half of the total material transport takes place by the time that D 
has decreased to 0^4 (at - 30°C), a mean Q (andY by Figure 18) is given by 


Q = Q n + A Q 


(7) 


£ Q T n - 30 
Qf -Q n " T n -T f 


( 8 ) 


since Q and T may both be assumed to vary linearly with time . With .f and ^ deter- 
mined F may be calculated by Equation (4), and ^determined by Equation (2) for 
assumed temperatures of nucleation . 

The nickel -profile in taenite approaches the < */X asymptotically, so the 
equilibrium concentration cai\not be conredlated with the kamacite composition . For 
practical purposes and for the resolution of the microprobe, Equations (2) and (4) were 
used to fix T^ at 350°C, corresponding to a C|| of 48.5 atom % nickel as plotted in 
Figures 2-U . 

Cooling model . The formula (Allan and Jacobs . 1956) for the central temperature T_ 
of a sphere of radius JR with thermal diffusivity at and surface temperature T cooling 
from an initial temperature T. is given by 
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T c ~ T s 
T i “ T s 


= 2 


Y2 ( ~ 


m + 1 - m Z ir 2 d t/ R 2 


1 ) 


m = 1 


( 9 ) 


2 2 

t r <x*t/R may be extracted with tables of I naersoll and Zobel (1948), for T n and 
and R calculated. T. was taken as I800°K, the minimum temperature at which liquid 
metal and silicate could have differentiated . c< was assumed to be 7 ♦ 10 cm /sec 
the approximate thermal conductivity for chondrites, as calculated from data compiled 
by Wood (1963). T £ =-IOO°K, the approximate temperature in the asteroid belt. 


- 19 - 


REFERENCES 

Agrell, S . O J . V . P . Long, and R . E . Ogilvie, Nickel content of kamacite near 
the interface with taenite in iron meteorites. Nature, 4882 , 749-750, 1963 . 

Allan, D. W. and J. A. Jacobs, The melting of asteroids and the origin of meteorites, 
Geochim . et Cosmochim . Acta, 9, 256-272, 1956. 

Allen, N . P . and C . C . Earley, The transformations <<.-*> V and in iron-rich binary 
iron-nickel alloys, J. Iron Steel Inst. (London), 166, 281-288, 1950. 

Anders, E., Origin, age, and composition of meteorites. Space Science Reviews, 
in press, 1964. 

Anders, E., and G. G. Goles, Theories on the origin of meteorites, J . Chem. Ed. , 

38, 58-66, 1961 . 

Armstrong, H. L., On solid-state diffusion with a linearly varying temperature, 

Trans. AIME , 212, 450-451 , 1958 . 

Bowles, J . S . and C . S . Barrett, Crystallography of transformations. Progress in Metal 
Physics, 3, 1-41, 1952. 

Carter, N . L. and G . C . Kennedy, Origin of diamonds in the Canyon Diablo and 
Novo Urei meteorites, J. Geophys . Res., 69 (12) , 2403-2421, 1964. 

Crank, J., The Mathematics of Diffusion, p. 146 , Clarenden Press, Oxford, 1956. 

Derge, G. and A. R. Kommel, The structures of meteoric irons, Am . J . Sci ., 34, 
203-214, 1937. 

Feller-Kniepme/er, M . and H . H . Uhlig, Nickel analyses of metallic meteorites by 
the electronprobe microanalyzer, Geochim. et Cosmochim. /^cta, 21, 257-265, 1961 . 

Fish, R. A., G. G. Goles, and E. Anders, The record in the meteorites, III . On the 
development of meteorites in asteroidal bodies, Astrophys . J ., 132 , 243-258, I960. 

Goldberg, E. , A. Uchiyama, and H. Brown, Tl.e distribution of nickel, cobalt, 


gallium, palladium and gold in iron meteorites, Geochim . et Cosmochim . Acta, 2, 
1-25, 1951 . 



- 20 - 


Goldstein, J . I The growth of the Widmanstatten pattern in metallic meteorites, 

Ph.D. thesis, Massachusetts Irstitute of Technology, 1964. Submitted to Geochim . 
et Cosmochim . Acta, 1964. 

Goldstein, J. I., R. E. Hannemand, and R. E. Ogilvie, Diffusion in the iron-nickel 
system at I atmosphere and at 40 kbar . , submitted to Trans . AIME, 1964 . 

Goles, G . G ., R > A . Fish, and E . Anders, The record in meteorites - I The former 
environment of stone meteorites as deduced from K -Ar ages, Geochim ■ et 
Cosmochim. Acta, 19, 177-195, I960. 

Hanson, D. and H. E. Hanson, The constitution of Ni-Fe alloys, J. Iron Steel Inst., 

102, 39-60, 1920. 

Henderson, E. P., El Burro, Coahuila, Mexico, meteorite. Am. Min., 26, 655-656, 

1941. 

Henderson, E . P . and S . H . Perry, A discussion of the densities of iron meteorites, 
G eochim. et Cosmochim. Acta, 6, 221-240, 1954. 

Henderson, E . P . and S . H . Perry, Studies of seven siderites, Proc . U . S . Natl . Museum, 
IQ 7 , 339-403, 1958. 

Hirano, K . M . Cohen, and B . L . Averbach, Diffusion of nickel into iron, Acta Met . , 

9, 440-4*5 , 1961 . 

Ingersoll, L . R., O . J . Zobel, and A . C . Ingerso!!, Heat Conduction, With Engineering 
and Geological Applications, Appendix H, McGraw-Hill, New York, 1948. 

Jones, F. W. and W. I . Pumphrey, Free energy and metastable states in the iron-nickel 
and iron-manganese systems, J. iron Steel Inst. (London), 163, 121-131, 1949. 

Jost, W., Diffusion in Solids, Liquids, Gases, pp. 68-75, Academic Press, New York, 
1952. 

Kaufman, L. and M. Cohen, The martensitic transformation in the iron-nickel system, 
Trans . AIME,* 206, 1393-1401, 1956. 

Kaufman, L. and M. Cohen, Thermodynamics and kinetics of martensitic transformations. 
Progress in Metal Physics, 7, 1 65-246, 1958. 



- 21 - 


Lovering, J. F., Structural and compositional studies on selected phases of iron and 
stony-iron meteorites, Ph, D. thesis, California Institute of Technology, 1956. 

Lovering, J. F., Differentiation in the iron-nickel core of a parent meteorite body, 
Geochim . et Cosmochim . Acta, 12 , 238-252, 1957a . 

Lovering, J. F., Pressures and temperatures within a typical parent meteorite body, 
Geochim. et Cosmochim. Acta, 12 , 253-261 , 1957b . 

Lovering, J. F ., The evolution of the meteorites — evidence for *-be co-existence of 
chondritie, achondritic, and iron meteorites in a typical parent m.’t«»cr*te body, 
in Researches on Meteorites , edited by Corleton B. Moore, pp. 179-198, John Wiley 
and Sons, New York, 1962. 

Lovering, J . F . and C . A . Andersen, to be submitted to Science , 1964. 

Lovering, J. F., W. Nichiporuk, A. Chodos, and H. Brown, The distribution of 

gallium, geranium, cobalt, chromium, and copper in iron and stony-iron meteorites 
in relation to nickel content and structure, Geochim. et Cosmochim. Acta, 2, 
263-278, 1957 . 

Maringer, R . E . , N . A . Richard, and A . E . Austin , Microbeam analysis of Widmanstatten 
structure in meteoric iron, Trans . AIME, 215 , 56-58, 1959. 

Mason, B., The origin of meteorites, J . Geophys . Res. 65 , 2965-2970, I960. 

Mason, & ., Meteorites , pp. 192-200, John Wiley and Sons, New York, 1962. 

Massalski, T . B ., Some metallurgical aspects of the study of meteorites, in Research es 
on Meteorites, edited by Corleton B . Moore, pp . 107-122, John Wiley and Sons, 

New York, 1962. 

Massalski, T . B . and F . R . Park, A quantitative metal lographic study of five octahedrite 
meteorites, J . Geophys . Res . , 67, 2925-2934, 1962. 

Meen, V . B ., The composition of the Maria Elena meteorite, Am . J . Sci ., 239, 412, 

1941 . 


- 22 - 


Meh 1 , R. F. and G. Derge, Studies upon the Widmanstatten structure, VIII the 
transformation in Fe-Ni alloys, Trans . AIMS, 125 , 482-500, 1937. 

Nichiporuk, W., Variations in the content o f nickel, gallium, germanium, cobalt, 
copper and chromium in the kamacite and taenite phases of iron meteorites., 

Geochim . et Cosmochim . Acta, 13, 233-247, 1958. 

Nichiporuk, W .-and H . Brown, The distribution of platinum and palladium metals in 

iron meteorites and in the metal phase of ordinary chondrites, submitted to J . Geophys . 
Res., July, i964. 

Owen, E . A., The structure of meteoritic iron, Phil . Mag ., 29, 553-567, 1940. 

Owen, E . A . and Y . H . Liu, Further X-ray study of the equilibrium diagram of the 
iron-nickel system,- J. Iron Steel Inst. (London), 163 , 132-137, 1949. 

Owen, E . A . and A . H . Sully, Equilibrium diagram of iron-nickel alloys, Phil . Mag. , 

27, 614-636, 1939. 

Patel, J . R. and M . Cohen, Criterion tor the action of applied stress in the martensitic 
transformation, Acta Met . , 1 , 53i-538, 1953. 

Perry, S . H ., The metallography of meteoric iron, U . S . Nat ■ Museum Bui I . , No . 184 , 
206pp. , 1944. 

Prior, G . T ., Catalogue of Meteorites, second edition, revised by M . H . Hey, 432 pp ., 
British Museum, London, 1953. 

Ringwood, A. E ., Chemical and genetic relationships among meteorites, Geochim. et 
Cosmochim. Acta, 24, 159-197, 1961 . 

Ringwood, A. E. and L. Kaufman, The influence of high pressure on transformation 
equilibria in iron meteorites, Geochim . et Cosmochim . Acta, 26, 999-1009, 1962. 

Uhlig, H ., Contribution of metallurgy to the origin of meteorites . Part 1-Structure of 


metallic meteqrites, their composition and the effect of pressure, Geochim. et 
Cosmochim. Acta, 6, 282-301, 1954. 




- 23 - 


Urey, H , C . , Diamonds, meteorites, and the origin of tiie solar system, Astrophys . J ., 124, 
623-637, 1956. 

Urey, H . C ., The early history of ihe solar system as indicated by the meteorites, 

Proc . Chem . Soc . (London) , 57-78, 1958. 

Urey, H . C ., Primary and secondary objects, J ■ Geophys . Res . , 64, 1721-1737, 1959. 

// 

Vogel, R., Uber die Struktur der Eiser-Nickel-Meteoriten, Z. anorg. u. allgem. Chemie. 
142 , 193-228, 1923. 

Vogel, R., Das System Eisen-Nickel -Phosphor und das meteorische Eisen, Handbuch der 
Metal Iphysik, II , 578-587, 1937. 

Vogel, R., Die Gefugeformen des Meteoreis ens und ihre Erklarung auf Grund des 
Zustandsdiagrammes des Systems Eisen-Nickel -Phosphor, Neues Jahrb . Mineral . 
AbhondL, 83 , 23-52, 1951 . 

Wood, J A., Physics and chemistry of meteorites. Chap . 12 in The Moon, Meteorites, 
and Comets, (Vol . !V of The Solar System), edited by B . Middiehurst and G . Kuiper, 

U . of Chicago Press, Chicago, 1963 . 

Yavnel* , A . A ., I . B . Borovski, N . P . Ilin, and I . D . Marchukova, A study of the 
composition of the phases of meteoritic iron by means of a local X-ray spectral 
analysis, Doklady Akad. Nauk S 3 S R, 123 , 256-258, 1958. 

Yeo, R . B . G . , The effects of some alloying elements on the transformation of Fe-22 .5% 

Ni alloys, Trans. AIME, 227 (4) , 884-890, 1963 . 

Young, J ., The crystal structure of meteoric iron as determined by X-rays, Proc . Roy . 

Soc. (London), 112 , 630-641, 1926. 



- 24 - 


TABLE I . Analysis of Standard Fe-Ni Alloys 
Weight % Ni Weight % Fe 


Chemical 

Probe 

Chemical 

(100-Ni) 

Probe 

5.15 

5.30 

94.85 

95.66 

10.26 

9.82 

89.74 

90.83 

15.10 

15.13 

84.90 

85.32 

19.25 

18.99 

80.75 

80.24 

25.17 

25.45 

74.83 

75.79 

36.46 

35 .62 

63.54 

65.52 

50.28 

49.07 

49.72 

51.31 



TABLE 2. Meteorites Studied 
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(5) Nichiporuk (1958) 
*Laue photograph diffuse 



TABLE 3 . Textural Classification of Iron Meteorites 
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LEGENDS FOR ILLUSTRATIONS 


Figure I . Photomicrograph of plessite band in Duchesne. Arrow shows typical position 
and direction of microprobe tracks . Scale bar 100^/ • 

Figure 2. Nickel concentration -distance profile across taenite band in Youndegin (Og). 

i 

— shows bulk nickel content . £ calculated amount of boundry movement , 

and x* are measured on section surface * 

Figure 3 . °rofile of nickel content across taenite band in Arispe (Og) . Peak nickel 
content in taenite established by diffusion calculation and numerous micro- 
probe analyses . 

Figure 4. Two profiles of nickel content in taenite and kamacite bands in Canyon Diablo 
(Og) . Narrow band shows plessite . 

Figure 5. Profiles of nickel across narrow and w T de bands in Toluca (Om) . 

Figure 6. Taenite band with central fine structure in Sacramento Mts . (Om). 

Figure 7. Effects of varying band-widths on taenite profiles in Spearman (Om). 

Figure 8. Profile near edge of fine plessite area in Duchesne (Of). Horizontal line 
shows average plessite concentration . 

figure 9. Profile across tine plessite band in Altonah (Of), showing average central 
plessite concentration . 

Figure 10. Profile across half of plessite area in Bristol (Of). 

Figure II . Variations of nickel and cobalt across a taenite band in Canyon Diablo having 
two coarse plessite lamcllas. 

Figure 12. Electron back-scatter picture at border of taenite and plessite (top) in Duchesne. 
Lighter areas are nickel-rich . Scale bar 

Figure 13. Electron back-scatter picture of nickel in interior of fine plessite in Duchesne. 
Lighter areas are taenite . Scale bar 5^ . 

Figure 14. Portion of Fe-Ni sub-solidus phase diagram at I atm . pressure, °C /KL + V" 
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boundry is uncertain . and T q refer to 2 ^ rans f ormc ^ on • 

Figure 15, Calculated cooling rate of meteorite during formation of Widmanstatten 
structure for different assumed kamacite nucleation temperatures (F n ) 
given by Figure 14. 

Figure 16. I - Cooling rate at center of sphere with initial temperature ISO^K, surface 

2 

temperature I00°K / end oC = 0.007 cm. /sec. 

2 - Effect of long-lived radioactivity . 

Figure 17. Calculated parent-body radius of meteorite at center during formation of 

Widmanstatten structure for different assumed kamacite nucleation temperatures 
(T ) given by Figure 14. 

Figure 18. Boundry-condition diffusion parameter as a function of bulk content C^ 
and boundry concentrations C ^ and in and phases . 
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